Astronomy & Astrophysics manuscript no. 13924 
May 27, 2010 



©ESO2010 



Adaptive optics near infrared integral field spectroscopy of 

NGC 2992* 

S. Friedrich^ R. I. Davies^ E. K. S. Hicks^ H. Engel\ F. Miiller-Sanchez^'^, R. Genzel^'^, and L. J. Tacconi^ 



' Max-Planck-Institut fiir extraterrestrische Physik, Postfach 1312, D-85741 Garching 

- Institute de Astrofisica de Canarias, E-38205 La Laguna (Tenerife), Spain 

^ Department of Physics, Le Conte Hall, University of California, Berkeley, CA, 94720, USA 



o 

(N 



(N 

o 
u 

6 

u 



in 

o 
o 



X 



ABSTRACT 

Aims. NGC 2992 is an intermediate Seyfert 1 galaxy showing outflows on kilo parsec scales which might be due either to AGN or 
starburst activity. We therefore aim at investigating its central region for a putative starburst in the past and its connection to the AGN 
and the outflows. 

Methods. Observations were performed with the adaptive optics near infrared integral field spectrograph SINFONI on the VUT, 
complemented by longslit observations with ISAAC on the VUT, as well as N- and Q-band data from the Spitzer archive. The spatial 
and spectral resolutions of the SINFONI data are 50 pc and 83 km s"' , respectively. The field of view of 3" x 3" corresponds to 450 pc 
X 450 pc. Bry equivalent width and line fluxes from PAHs were compared to stellar population models to constrain the age of the 
putative recent star formation. A simple geometric model of two mutually inclined disks and an additional cone to describe an outflow 
was developed to explain the observed complex velocity field in H2 l-OS(l). 

Results. The morphologies of the Bry and the stellar continuum are different suggesting that at least part of the Bry emission comes 
from the AGN. This is confirmed by PAH emission lines at 6.2 yum and 11.2/im and the strength of the silicon absorption feature 
at 9.7 /im, which point to dominant AGN activity with a relatively minor starburst contribution. We find a starburst age of 40Myr 
- 50 Myr from Bry line diagnostics and the radio continuum; ongoing star formation can be excluded. Both the energetics and the 
timescales indicate that the outflows are driven by the AGN rather than the starburst. The complex velocity field observed in H2 1 -0S( 1 ) 
in the central 450 pc can be explained by the superposition of the galaxy rotation and an outflow. 

Key words, galaxies: active - galaxies: individual: NGC 2992 - galaxies: Seyfert - galaxies: starburst - infrared: galaxies 



1. Introduction 

NGC 2992, initially classified as a Seyfert 2, was later classified 
as an intermediate Seyfert 1 galaxy on the basis of a broad Ha 
component with no corresponding li/3 component in its nuclear 
spectrum (Ward et al. |1980l ). Later, Glass ( 119971 ) suggested that 
it is a hybrid between an intermediate Seyfert and a starburst 
galaxy, induced by the interaction with NGC 2993 to the south. 
Both galaxies are connected by a tidal tail with a projected length 
of 2:9. 

NGC 2992 is highly incUned by about 70° to our Une 
of sight with a broad disturbed lane of dust in its equatorial 
plane. The redshift is 0.0077 which corresponds to a distance 
of 32.5 Mpc and an angular scale of 150pc per arcsec assuming 
i/o=75kms-'Mpc-'. 

At a wavelength of 6 cm NGC 2992 shows prominent loops 
to the northwest and southeast. This feature, designated as 
figure-of-8, has a size of about 8" with the long axis at a po- 
sition angle (PA) of -26°, which is not coincident or perpen- 
dicular to the major axis (Ulvestad & Wilson 1984). Chapman 
et al. (|2P00) conclude from their IR observations that the most 
convincing interpretation for these loops is that of expanding gas 
bubbles which are seen as limb-brigthened loops as was already 
proposed by Wehrle & Morris ( 119881 ). Two possible explanations 
are discussed for the origin of the loops: They may be driven by 



an AGN or a superwind from a compact starburst. X-ray ob- 
servations point to the first possibility: Colbert et al. ( 1998 ) con- 
cluded that the large scale outflows are driven by nonthermal jets 
from the AGN that entrain the suiTounding material and heat it 
by shocks over kpc spatial scales. 
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at the European Southern 



Heckman et al. (1981') reported complex kinematics in NGC 
2992 on the basis of longslit spectra at R4= 120°. This work was 
extended by Colina et al. (1 19871 ) who find indications for an out- 
flow within a plane highly inclined to the plane of the galaxy. 
Marquez et al. (119981 ) measured [Oiii] and Ha lines at nine 
different position angles and concluded, although from a sim- 
ple kinematical model, that the double peaked line profiles and 
asymmetries can be accounted for by an outflow within a conical 
envelope or on the surface of a hollow cone. The high spectral 
resolution observations (4.9 A: 3440 A - 5579 A, 1 1 .9 A: 5840 A 
- 9725 A) of naiTow line regions in NGC 2992 of Allen et al. 
(1999) show at most positions a double component line profile. 
One component of the double peaked lines follows the galactic 
rotation curve, while the other is identified as a wind that is ex- 
panding out of the plane of the galaxy with a velocity of up to 
200 kms'. These authors did not find spectral evidence for a 
significant contribution from a starburst to the outflow. 

Depending on the inclinations of galaxy and cone several 
interpretations of these kinematics are possible. It is generally 
assumed that the NW edge of the disk of NGC 2992 is closest 
to us with trailing spiral arms, and that the SE cone is closer 
and directed at us. However, an arc-shaped emission found by 
Garcia-Lorenzo et al. (120011 ). which might be associated with 
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the radio loop, is red shifted in the SE, and can therefore only be 
consistent with an outflow if it is pointing away from us. 

In this paper, we present spatially and spectrally resolved K- 
band observations of the inner 3" x 3" of NGC 2992 taken with 
the adaptive optics integral field spectrograph SINFONI. These 
are complemented by longslit observations along the major axis 
out to a radius of about 5" with ISAAC, and N- and Q-band 
spectra from the Spitzer archive. SINFONI data were taken with 
unprecedent spatial resolution in the K-band which allows us to 
study the dynamics of gas and stars in the central region in de- 
tail. Together with the spectral resolution of SINFONI this en- 
ables us to search for a starburst and its connection to the AGN 
and, possibly, to the outflows. After presenting the observational 
results in section 3, we discuss the observations in section 4, and 
show that there is evidence for a nuclear starburst in NGC 2992. 
Finally, in section 5 we concentrate on the dynamics of the gas 
and stars and the consequences on the geometry of the galaxy. 



2. Observations and data reduction 

2.1. SINFONI 

K-band spectra of NGC 2992 between approximately 1 .95 fim 
and 2.45 ;um were taken on 13 March 2005 at the VLT with 
SINFONI, an adaptive optics near infrared integral field spectro- 
graph (Eisenhauer et al.'2003, Bonnet et al. 2004") with a spectral 
resolution of 3600. The pixel scale was 0'.'05 x (J.' 1 and the ex- 
posure time 600 s for both object and sky frames. Object frames 
were interspersed with sky frames using the sequence O-S-O-O- 
S-O. A total of 12 object frames were obtained, resulting in a 
total exposure time of two hours. For a near diffraction limited 
correction the nucleus of the galaxy was used by the AO module. 

Data were reduced in a standard manner by subtraction sky 
frames, flatfielding, correcting bad pixels, and wavelength cal- 
ibration with an arc lamp using the SPRED software package 
(Abuter et al. 2006), which also performs the reconstruction of 
the data cube. No spatial smoothing by a median filter was ap- 
plied. Telluric correction and flux calibration were performed 
with the G2V star HD 85724 (K=7,881). Residuals from the 
OH line emission were removed using the method described in 
Davies (I2007al l. 

The spatial resolution has been estimated from both the 
broad Bry and the non-stellar continuum (Davies l2007bi Both 
methods yield symmetric PSFs with FWHMs of a.'32 and a.'29, 
respectively, corresponding to 48 pc and 44 pc, respectively. The 
spectral resolution derived from sky lines is 83kms '. NGC 
2992 is part of a (heterogeneous) sample of type 1 and type 2 
Seyfert galaxies, ULIRGs, and a Quasar, which were observed 
to investigate the connection between star formation and AGN. 
Further details on data reduction can be found in Davies et al. 
(|2007] | and Hicks et al. (|2009] |. 

2.2. ISAAC 

NGC 2992 was also observed on 13 January 2006 with the 
Infrared Spectrometer And Array Camera ISAAC at the VLT 
in order to obtain wide FOV longslit spectra centered on 2. 1 fim 
(to measure H2 l-OS(l)) and 2.29 yum (to measure '^00(2-0)). 
The 2' long and 0'.'3 wide slit was aligned parallel to the 
morphological large scale major (PA=34.62°) and minor axes 
(PA=-55.38°), respectively. The exposure time amounts to 40 
minutes per slit position. The spectra were reduced using stan- 
dard techniques carried out by the ESO pipeline software. 
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Fig. 1. SINFONI-Specti-um of NGC 2992 integrated over the 
whole FOV of 3" x 3". The most prominent lines are those 
of H2l-0S(3) at 1.958//m, [Sivi] at 1.96;um, H2 l-OS(l) at 
2.12/im and Bry at 2.16//m. 

Table 1. Measured line fluxes for NGC 2992 in a field of view 
of r.'5 X r.'5 centered on the nucleus 



Line 



r Fiux^ Flux;. 

(pm) (10-"*Wm-2) (10-'**Wm--) 



H2 1-0S(3) 


1.9576 


11 ±0.2 


5 


[Si VI] 


1.9634 


13 ±1.6 


4 


H2 1-0S(2) 


2.0338 


3.9 ±0.7 


24 


H2 l-OS(l) 


2.1218 


9.4 ±0.2 


2.5 


Bry 


2.1661 


14.3 ±3.5 


11.6 


Bry (narrow) 


2.1661 


4.8 ±1 


- 


H2 l-OS(O) 


2.2233 


24 ±1.2 


- 



" wavelengths are given in the rest frame 

* this paper 

'■ from Gilli et al. <2000b . aperture size 1 " x 2 " 



In order to determine line positions, and subsequently ra- 
dial velocities of H2 l-OS(l), a spectrum from the outer edge of 
the galaxy's major or minor axes was extracted, the signal-to- 
noise determined and if necesary the spectrum from the adjoin- 
ing pixel added until a S/N=3 was reached. This procedure was 
repeated along the slit to the opposite edge of the galaxy result- 
ing in a set of 1 14 and 62 individual spectra along the major and 
minor axes, respectively. Bry (2.16//m) can also be detected in 
these spectra. These 2D-data sets were then analysed in a similar 
way to their 3D counterparts (see section 5). 

3. Observational results 

3.1. The spectrum of NGC 2992 

In Fig. [U we show the SINFONI spectrum of NGC 2992 inte- 
grated over the whole FOV of 3" x 3" centered on the nucleus. 
The most prominent lines are H2 1-0S(3) at 1.958jum, [Sivi] 
at 1.96yL/m, H2 l-OS(l) at 2.l2fj.m, and Bry at 2.16jum, other 
molecular hydrogen transitions and some CO bandheads are also 
present. The line fluxes are summarized in Table [1] together with 
Une fluxes from GiUi et al. (2000). 

Variability of the infrared flux of NGC 2992 was already re- 
ported by Glass ( 1997). Apart from an outburst in 1988, which 
was especially strong in the K- and L-band, he reported a fading 
of NGC 2992 by about 20% in the K-band, 50% in the L-band 
between 1978 and 1996, and by a factor of w 20 in X-rays at the 
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Fig. 2. Continuum, velocity, and dispersion maps of NGC 2992 (left to right) of stars, Bry, and H2 l-OS(l) (top to bottom). 1" 
corresponds to 150 pc. The contour lines of the Bry continuum are overlayed on all continuum maps and a cross indicates the center 
of each map. North is up and east to the left. 



same time. Then, in 1998 BeppoSAX observations showed the 
X-ray flux to be again similar to that of 1978 (Gilli et al. 2000). 
Also line properties of Balmer- and Paschen-lines in the near 
infrared were similar to those observed before the decline. It is 
obvious from Table[T]that our line fluxes from 2005 (determined 
for a similar aperture size) are higher than those of Gilli et al. 
(I2OOOI 1 from January 1999. Unfortunately we have no informa- 
tion about Ha and Pa/3 to compare them to former observations 
given in Gilli et al. 

Oliva et al. ( fT995l l show that the '^00(2-0) at 2.2//m in 
the K-band, and the '^CO(6-3) bandhead at 1.6yum in the H- 
band, can be used to separate stellar and non-stellar continua 
in AGN. However, doing so requires high signal-to-noise data 
in both bands. Instead, we follow the method of Davies et al. 
(I2007I I which uses just the '^CO(2-0) bandhead to calculate the 
fraction of the nuclear flux that is stellar By assuming that the 
stellar equivalent width in the nuclear region is about 12 A, re- 
gardless of the star formation history, as predicted by the mod- 
els created with the stellar synthesis code STARS (Sternberg 
11998 ; Sternberg et al.^003), which is comparable to Starburst99 
(Leitherer et al. 1 19991 Vazquez & Leitherer l2005l) . we calculated 
the non-stellar dilution using the equation: 



1 -D 



£Wi. 



where D is the fraction of the continuum that is non-stellar We 
found an equivalent width of (5.0+ 1 .6) A in the central 0'.'5 xO'.'5, 



implying that 58% of the K-band luminosity in this region comes 
from non-stellar components. With this value we determined a 
stellar luminosity of Lk=3.9-10^ Lq and a total continuum lumi- 
nosity of L™'"=5.4-10^ L0 in this central region in the K-band. 



3.2. 2-climensional flux distributions, velocities, and velocity 
dispersions 

2-dimensional flux distributions, velocities, and velocity disper- 
sions of individual lines were determined by fitting the unre- 
solved line profile of an OH line, which was convolved with 
a Gaussian, to the respective lines. The continuum level was 
found by a linear fit to the spectrum and subtracted. The uncer- 
tainty of a fit was estimated using x^ techniques by refitting the 
best-fit Gaussian with added noise at the same level as the data 
100 times. Standard deviations of the best-fit parameters are then 
used as the uncertainties. Uncertainties of velocities and veloc- 
ity dispersions are of the order of 10 km s"' to 15 km s"' for the 
'^CO(2-0) bandhead at 2.2/im (i.e. for the stars) and Bry, and 
about 5kms"' for H2 l-OS(l). Errors of the continuum range 
from below 1% for Bry and H2 l-OS(l), and 6% for the stars. 

In Fig.|2]we show the morphology of the stellar continuum, 
the Bry- and H2 l-OS(l) lines and their respective velocities and 
velocity dispersions. The spatial resolution is 0.3"corresponding 
to a length scale of about 50 pc; the spectral resolution is 
83kms '. The continuum seems to trace an inclined disk with 
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fey- ^'isr'i 



Fig. 3. Velocity (left) and velocity dispersion maps (right) of H2 l-OS(l) with the contours of the radio flux at 8.46 GHz (Leipski et 
al.'2006' courtesy Heino Falcke) and the rough positions of the redshifted and blueshifted velocity residuals of [Oiii] (dashed lines 
at bottom and top) found by Garcia-Lorenzo et al. (2001) overlayed. North is up and east to the left. 
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Fig. 4. Supernova rate and Bry equivalent width calculated with STARS as a function of age for several star formation timescales. 
They are normalised to the K-band stellar continuum and Lk is the total luminosity in the 1.9-2. 5 jum band in units of the solar 
luminosity (Lio/=3.8xlO^^W). The observed values are drawn as horizontal lines, the error margin is indicated by dashed lines. 



the north-west side more obscured. Given the broad dust lane 
seen in optical images this is not surprising. 

There exist remarkable differences. Firstly, the morphology 
of the Bry line does not follow the stellar continuum. To make 
this more obvious the contour lines of the Bry continuum are 
overlaid in Fig. |2] on the stellar continuum. The center of the 
stellar continuum is shifted to the south east and the intensity 
distribution is elongated in the north east to south west direction 
in contrast to the circular Bry contour. This means that part of 
the Bry emission must be attributed to other sources than a star- 
burst. On the other side, the morphology of the H2 l-OS(l) line 
is similar to the stellar continuum. 

The zero line of the velocity is orientated roughly east west 
in all three velocity maps. The maximum velocities range from 
H-lOOkms"' to -nlTOkms"' and the minimum velocities from 
-100 km s"' to -140kms"'. The velocity map of the stars shows 



a distinct increase to up to 100 km s 
of the Bry continuum. 



at the position of the centre 



The stellar component, Bry and H2 l-OS(l) all show veloci- 
ties away from us in the southwest and towards us in the north- 
east. The velocity field of H2 l-OS(l) and also that of Bry, albeit 
not as clearly, show redshifted emission also to the southeast. 
Garcia-Lorenzo et al. (,2001 ) investigated the velocity field of 
[O III] and found an arc of redshifted velocities to the southeast 
and one of blueshifted velocities to the northwest with velocities 
of (64 + 18)kms"', and (-79 + 28)kms"', respectively. These 
arcs are just outside our field of view. Only the redshifted ve- 
locities in the south east corner of H2 might be caused by them 
(Fig-EJ. Since these "arcs" resemble the loops of the figure-of-8 
when superimposed on the radio emission there might be a re- 
lation between both structures. However, this would suggest an 
inflow, if the southeast part of the figure-of-8 is closer to us, as 
generally assumed. In the picture of expanding gas bubbles this 
cannot be understood. 

Already Chapman et al. (2000) proposed that the radio mor- 
phology consists of a component out of the plane of the galaxy 
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disk and one within the disk associated with a spiral arm. On 
the basis of the narrow and well defined spectral lines within the 
arcs Garcia-Lorenzo et al. conclude that the arcs are due to gas 
in the disk rather than in the loops. However, none of the arcs 
is associated with the structures that Chapman et al. interpret as 
spiral arms. 

There is a small region in the velocity field of Bry which ex- 
tends from west of the centre to nearly the edge of the FOV in 
the SW with a lower velocity relative to its surroundings. This 
can be interpreted as an outflow superimposed on the galaxy ro- 
tation. Although it coincides with the southern inner spiral arm 
found by Chapman et al. (2000) in their model-subtracted R- 
and H-band images, it cannot be related to it because then the 
spiral arm would move relative to the ambient medium. There is 
no counterpart in the H2 l-OS(l) and stellar continuum velocity 
fields. 

Common to both the stars and the lines are relatively high 
velocity dispersions of up to 250 km s"' for the stars and Bry, 
and 130km s~^ for H2 l-OS(l). However, the stellar dispersion 
map shows no structures, except that there seems to be a local 
minimum at the position of the centre of the Bry continumm. 
There, the measured velocities range from (110 - 150) km s"' 
±(20 - 30)kms"' compared to the surrounding values of more 
than (200 + 15)kms '. This high stellar velocity dispersion 
which we observe even on the smallest scales suggests that the 
K-band light is dominated by stars in the bulge. 

The high velocity dispersion of Bry is confined to the inner 
1'.'2, outside dropping to (0 - 30) km s ' . The dispersion of H2 1 - 
0S(1) shows maxima of 150 km s"' to the southwest and north- 
east of the centre, which correspond with the velocity maxima 
and minima. They also lie at the edge of radio features within 
the loops (Fig. [3]) and may be associated with inner spiral arms. 
The dispersion of H2 l-OS(l) also shows a local minimum at the 
centre where the dispersion drops to 70 km s ' . There might also 
be a local maximum at the position of the centre of the Bry con- 
tour 

The complexity and the differences of velocity and disper- 
sion maps of the stars and the lines might be due to the super- 
position of the galaxy rotation and an outflow as it was already 
suggested by several authors (e.g. Chapman et al. 120001) . 

4. Constraining nuclear star formation in NGC 2992 

The possible existence of a nuclear starburst in NGC 2992 was 
investigated by Davies et al. ( 120071) . who compared the radial 
profile of the stellar continuum to an r^^'^ law and an exponen- 
tial profile. They found that there might be excess continuum 
at r < 0'.'5 due to a distinct stellar population superimposed 
on the bulge. However, although the size scale matched the nu- 
clear star forming regions in other AGN with typical HWHM of 
10pc-50pc (s. Fig. 5 in Davies et al.'2007), they did not reach a 
conclusion on NGC 2992 because the excess was only apparent 
for an exponential profile. Data at a higher resolution of about 
(X'l would be needed to explore the spatial properties further 
Instead, here we look at various star formation diagnostics in 
order to assess the luminosity and age of a putative nuclear star- 
burst. 

The stellar population synthesis code STARS (e.g. Sternberg 



11998 i Sternberg et al. ,2003J was also used to model the equiv- 
alent width of Bry and the supernova rate. The program fol- 
lows the evolution of a cluster of stars through the Hertzsprung- 
Russell-Diagram depending on the star formation history. We 
assumed decaying star formation rates with timescales of 10*' 
(instantaneous), 10^, 10^, and 10'^ (continuous) years. In Fig. 



|4] the results from STARS together with the observations are 
presented. Further details on star formation diagnostics can be 
found in Davies et al. (i2007j . 

4.1. Bry equivalent width and supernova rate 

Using the stellar continuum luminosity from section 3.1, an up- 
per limit for the equivalent width of Bry associated with star for- 
mation from the narrow Bry line flux can be estimated. However, 
for NGC 2992 it is difficult to quantify what fraction of Bry is 
associated with star formation, since the morphology of the line 
does not follow the stars (Fig.|2]i. Especially the south-west side 
shows velocities which are blue shifted relative to their environ- 
ment indicative of motions towards us. The different morphol- 
ogy of H2 l-OS(l) and Bry found in section 3.2 points towards a 
kinematically different origin, and supports the conclusion that 
at least part of the Bry emission is due to the AGN and not to star 
formation. Thus the narrow Bry line flux includes a certain un- 
known fraction from the AGN and the derived equivalent width 
of (7.5 + 0.5) A is only an upper limit. This value, and all the 
numbers throughout the rest of the paper, are determined in the 
central 0'.'5 x 0'.'5 unless otherwise specified. 

We estimate a supernova rate for NGC 2992 from the un- 
resolved radio flux of 7 mJy at 5 GHz (Wehrle & Moms [19881 ) 
and the estimation done by Sadler et al. (1 19951 ) that the core flux 
at 5 GHz is less than 6 mJy. Their estimation is based on their 
higher spatial resolution measurements at 2.3 GHz and the non- 
detections at 1.7 GHz and 8.4GHz. Taking a flat spectral index 
(Chapman et al. 120001 ) the 5 GHz flux will not be much less than 
6 mJy, which would result in about 1 mJy extended emission. 
If we assume that it can be attributed to star formation, a su- 
pernova rate of 0.003 yr ' can be derived with the relation of 
Condon d 19921 ). On the other hand, if we assume a spectral in- 
dex of 0.82, which was derived by Ulvestad & Wilson ( 119841 ) 
for the unresolved flux at 1 .5 GHz and 5 GHz, together with the 
2.3 GHz flux of Sadler et al. (T995 ) we get a flux of about 3 mJy 
at 5 GHz. This leaves us with about 4 mJy extended emission 
corresponding to a supernova rate of 0.011 yr '. We consider 
the derived supernova rates as extremes and adopt a mean value 
of 0.007 yr-i. 

We can now try to determine a star formation age from 
the observed Bry equivalent width and the supernova rate by 
comparing them to the theoretical values calculated by STARS, 
which are plotted in Fig.|4]for different star formation timescales. 
It is immediately obvious that for a continous star formation 
the observed Bry equivalent width and the supernova rate yield 
inconsistent star formation ages of more than 3 Gyr and (25 
- 50) Myr, respectively. In addition, an ongoing star formation 
would require an equivalent width of Bry of 10 A - 15 A, which 
is not observed. 

An intermediate decay timescale of 10** yr yields a star for- 
mation age from the Bry equivalent width of about 280 Myr and 
25 Myr - 60 Myr from the supernova rate, which is an order of 
magnitude lower. Even if we consider the lower boundary of the 
supernova rate of 0.003 yr' the resulting star formation age of 
160 Myr is well below the age derived from the Bry equivalent 
width. 

The best agreement between the derived star formation ages 
from the Bry equivalent width and the supernova rate can be 
achieved for a decay timescale of 10^ yr. The resulting star for- 
mation ages amount to 40 Myr and 50 Myr respectively. The 
largest error of + 10 Myr is due to the uncertain radio flux at 
5 GHz and subsequently the supernova rate. The unknown con- 
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tribution of the AGN to Bry does not contribute significantly to 
the uncertainty, because the curve is steep. Even an increasing 
contribution of the AGN up to nearly 100%, which leads to a 
respective lower stellar contribution to the Bry equivalent width, 
would increase the age only up to 60 Myr. This is still within the 
errors set by the supernova rate. We therefore conclude that a 
stellar population exists in NGC 2992 with an age of 40 Myr to 
50 Myr. 

In principal there is also agreement between the Bry equiv- 
alent width and the supernova rate for instantaneous star forma- 
tion at a very young age of about 10 Myr, but we consider this 
less likely; partly because the timescales are much shorter, but 
primarily because it requires very careful tuning to match the in- 
creasing supernova rate and the quickly decaying Bry equivalent 
width. 

Combining the age range of 40 Myr - 50 Myr with the bolo- 
metric luminosity and black hole mass for NGC 2992 allows us 
to locate this galaxy on Figure 11 of Davies et al. (2007). This 
figure shows how the luminosity of an AGN might be related 
to the age of the starburst, and NGC 2992 can be placed in the 
lower left where the starburst age is young and the AGN is, at 
the present time, accreting at a rate that is low compared to the 
Eddington rate. The galaxy thus appears to fit into the proposed 
scheme, where there is either a delay between the onset of star- 
burst activity and the onset of AGN activity, or star formation 
ceased once the black hole become active (s. Davies et al. 120071 
for a detailed discussion). 

We also briefly consider whether the interaction between 
NGC 2992 and NGC 2993 might be responsible for the star- 
burst. The gravitational (tidal) force during an interaction can be 
approximated as a delta-function that occurs at perigalacticon. 
According to Due et al. (I2000I I perigalacticon of NGC 2992 and 
NGC 2993 occurred about 100 Myr ago. At that time the gas 
would have felt the strongest impulse. Given that it takes time 
for the gas to respond and, for example, be driven to the centre, 
it is perhaps plausible for the interaction to have triggered the 
starburst. 



4.2. Mass-to-light-ratio and bolometric luminosity 

In order to estimate the mass-to-light-ratio we determined the 
dynamical mass within 0'.'5 from the following relation; 



Mdyn = (vi, + 3(r2)-/?/G 



(2) 



where Vrot is the observed velocity corrected by the inclination 
and cr is the observed velocity dispersion at radius R (Bender et 
al. |1992t Davies et al. 2007). We did not correct for the molecu- 
lar gass mass nor for the contribution of an old stellar population. 
Therefore the derived dynamical mass of 8T0^ Mq is only an 
upper limit. With a black hole mass of 5-10^ Mq determined 
from velocity dispersion measurements (Woo & Urrv l2002l l the 
dynamical mass is at most 3-10^ Mq for the young popula- 
tion. Taken at face value, the resulting mass-to-light-ratio of 0.8 
speaks in favour of a star formation age greater than 80 Myr for 
all star formation rates. This is a factor of 2 higher than the age 
derived from the Bry equivalent width and the supernova rate. 
However, it needs to be borne in mind that the dynamical mass is 
only an upper limit of the younger stellar population (it will con- 
tain a significant contribution from the older bulge population) 
and the uncertainty of the black hole mass can be nearly as high 
as 50% (s. Woo & Urry 2002 ). Either a change of the black hole 
mass by 10% or a reduction of the mass of the young population 
by 20% would bring the star formation ages into agreement. 



STARS also allows us to estimate the stellar bolometric lu- 
minosity. The ratio between Lboi and L^ depends on the age and 
the exponential decay timescale of the star formation. Given the 
star formation age of 40 Myr - 50 Myr we derive a bolometric 
luminosity for NGC 2992 of 70-100 times the K-band lumi- 
nosity or (3-4)-10^ L0, which is (6-8)% of the total bolometric 
luminosity of NGC 2992 (5-10'° L© calculated between 8 and 
1 000 fim from IRAS 1 2 - 1 00 fim flux densities, and additionally 
taking into account the emission at shorter wavelengths; Davies 
etal. l2007l l. 

4.3. Star formation rate 

Using the star formation history (decay timescale and age) 
above, we can estimate the star formation rate. The first step is to 
determine what fraction of the K-band luminosity in the central 
0'.'5 comes from young stars. To do so, we estimate the contri- 
bution from an old population, by assuming that the dynamical 
mass is dominated by this old population. For an adopted age of 

5 Gyr, STARS yields a mass-to-Hght ratio of about 3OM0/La:. 
The dynamical mass determined above of 3 ■ 10^ Mq then results 
in a luminosity of 1 ■ 10^ Lq for the old population, one order 
of magnitude below the observed K-band luminosity. Therefore 
the K-band luminosity is dominated by the young population. 
We can therefore scale the luminosity of the starburst model to 
the total K-band luminosity in the central 0'.'5, yielding an initial 
star formation rate of 4. 3 MQyr"' for an age of 50 Myr. Since the 
SFR is decaying, an alternative, perhaps more meaningful, num- 
ber may be the time-averaged star formation rate - which we de- 
fine as simply the mass of stars formed divided by the age. This 
yields (SFR) ~ 1 MQyr"'. While this appears modest, in terms 
of SFR per unit area it corresponds to ~ 200MQyr"'kpc"^, 
within the range of values found for the nuclear starbursts in 
other nearby AGN by Davies et al. ( I2007I I. 

An independent method to estimate the SFR is via the poly- 
cyclic aromatic hydrocarbon (R4H) features. Farrah et al. (120071) 
employed diagnostics based on the luminosities and equivalent 
widths of fine-structure emission lines and R4H features as well 
as the strength of the 9.7 jum silicate absorption feature to in- 
vestigate the power source of the infrared emission in ULIRGs. 
With these tools they were able to distinguish between starburst 
and AGN caused infrared emission, and to determine a star for- 
mation rate. 

In order to apply these diagnostics to NGC 2992 we used 
near infrared pipeline reduced IRS spectra from the Spitzer 
archive. If more than one spectrum was available we took the 
mean value from all spectra to determine line luminosities. 
25 i^m and 60 /im fluxes were taken from the IRAS Revised 
Bright Galaxy Sample (Surace et al . 120041 1. The results are given 
in Table |2] 

The star formation rate 

SFR [MQyr-'] = 1.18 x 10-4'(i6.2;,m + iii.2;,m) [erg s"'] (3) 

was empirically derived from the PAH 6.2 yum + 11.2/im lumi- 
nosity for a star formation occurring on timescales of 10^ - 10** 
years (Farah et al. 120071 ). It yields 2.5 Mq per year or about 
lOMQyr 'kpc"2, if we take the slit width of 3'.'6 of the IRS 
spectra as the size of the investigated region for NGC 2992. This 
value is consistent with the SFR derived above, being above the 
time averaged rate but below the initial peak rate. While it can be 
regarded as fairly normal compared to the star formation rates of 
100 Mq per year for nearby luminous starbursts (e.g. Weedman 

6 Houck 120081 ). in terms of SFR per unit area, it is as extreme 
as these luminous starbursts. 
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Table 2. Measured line luminosities of fine-structure lines and 
PAHs for NGC 2992 



injected at a constant rate and expanding into a uniform medium 
with a density no in cm"^ (Heckman et al. 119901 1. From 



Line A" luminosity 
(fim) (10^3 W) 



dE/dt « 3 X 10*'r2p^v^oo«o ergs s"' 



(4) 



PAH 


6.2 


9.2+0.7 


PAH 


11.2 


12.5±0.2 


[Siv] 


10.5 


2.0+0.01 


[Siii] 


18.7 


3.8±0.3 


[Neil] 


12.8 


6.8±0.9 


[Nelll] 


15.6 


7.1±0.4 


[Nev] 


14.3 


3.9+1.1 


[Oiv] 


25.9 


14.7+1.6 



" wavelengths are in the rest frame 

4.4. Starformation versus AGN contribution 

The strength of the 7.7 ;um PAH feature (1.7+0.02) indicates that 
the AGN in NGC 2992 contributes more than 50% to the infrared 
flux between 1 ;um and 1000 /im (Genzel et al. 1998). [Ne V] at 
14.32 yum, which is strong in AGN spectra but not in star form- 
ing regions, and [Oiv], which is strong in star forming regions, 
are both present in the spectra of NGC 2992, pointing to a co- 
existing starburst and AGN activity. 

In order to quantify this a little more we followed Farrah 
et al. (120071 1 and compared the equivalent width of the PAH at 
6.2//m(EW=(0.084 + 0.002);um)totheratioof [Nev] to [Neii] 
and [Oiv] to [Neii], respectively. From these diagnostics we de- 
rive a 10% - 20% starburst fraction and about 50% AGN fraction 
to the infrared flux for NGC 2992 in agreement with the result 
from the PAH 7.7 //m feature. The AGN contribution can also 
be estimated from the ratio of the 25 yum flux to the 60 /im flux 
against the [Ne v] to [Ne ii] ratio. Again the AGN contribution 
amounts to about 50%. 

In this and the previous section we utilized different methods 
to estimate the contribution of stellar, non-stellar and starburst 
components to the K-band and infrared flux (1 fim - IQOQfim). 
Independent of wavelength range and FOV we found a 50% - 
60% AGN contribution. From the CO bandhead at 2.2 fim we 
found that 42% of the K-band flux in the central 0'.'5 xO'.'5 can be 
attributed to stars and from the star formation history that these 
stars must be young. Also in the 3" x 3" FOV the K-band light is 
dominated by stars. However, according to the high stellar veloc- 
ity dispersion it is more likely that stars from the bulge dominate 
in this larger FOV. Finally, we find from PAH line diagnostics in 
the infrared a 10% - 20% contribution to the infrared flux from 
a starburst. In addition to the AGN and starburst components, 
there still remains about 20% - 40%, which must be attributed 
to a third component - most likely the disk and bulge of the 
galaxy. This means that while the young starbust is important in 
the central tens of parsecs, it is still only a small part of the in- 
tegrated luminosity of the galaxy. Therefore we conclude that in 
NGC 2992 we see a dominating AGN activity with a relatively 
minor contribution from a starburst. 



5. Origin of nuclear outflow 

5.1. Starburst or AGN driven outflow 



Heckman et al. (1 19901 1 provided evidence that galaxies with a 
central starburst have large scale mass outflows which are pre- 
sumably driven by the kinetic energy supplied by supernovae 
and winds from massive stars. In order to examine whether the 
outflow in NGC 2992 can be caused by the starburst we calcu- 
lated the radius r in kpc of a bubble which is inflated by energy 



wifli the velocity (in lOOkm s"') v iop w 2 (Allen et al. 1 19991 1. 
no ~ 3 cm"-^ (Heckman et al. '1990\ and dE/dt equal to the me- 
chanical luminosity from stellar winds and SNe taken from star- 
burst99 for instantaneous (10'*" ergs"') star formation we get a 
radius of 40 pc. This radius is too small compared to the ob- 
served radii of about 600 pc of the figure-of-8 and 2 kpc for the 
outflows found in He and [Oiii] (Allen et al. ll999l l. An energy 
injection of 3 • lO'*^ erg s ' would be required to explain the ob- 
served radius and velocity of the outflow in NGC 2992. This 
estimation is based on a symmetric geometry. If we instead as- 
sume a cone with an opening angle of 1 16°, as seen for the large 
scale outflow (Allen et al. 1999), the energy needed is a factor of 
10 lower but still cannot be provided by the starburst. 

The 2 kpc extension of the large scale outflow and the veloc- 
ity of 200kms"' suggest that the outflow originated in an event 
at most 10 Myr ago. This is inconsistant with our derived star 
formation age of (40 - 50) Myr but in agreement with the typical 
timescale of the active phase of a black hole in an AGN. We note 
that although it is also consistent with the age for instantaneous 
star formation, this scenario we regarded as rather unlikely since 
it requires fine-tuning (see section 4.1). 

Thompson et al. (2005) showed that optically thick pres- 
sure supported starburst disks have a characteristic flux of about 
lO'^'Lgkpc"^. Above this limit matter would be blown away. 
With the derived luminosity of (3 - 4) • 10^ Lq within the cen- 
tral 0'.'5 we estimated the characteristic flux for NGC 2992 to 
be (5 - 7) ■ 10" LQkpc"^ which is well below this limit. Only 
when the starburst was at its peak luminosity, which for our pre- 
ferred model is nearly 10 times greater than the current lumi- 
nosity, would the star formation approach - although remaining 
below - this threshold. The conclusion from both energy and 
timescale considerations, is that for the starburst to drive the out- 
flows it would have to be more intense and more recent - both 
of which are ruled out by our analysis. 

Instead, we consider whether the outflows in NGC 2992 can 
be driven by the AGN. In order to estimate the AGN luminos- 
ity of NGC 2992 we use the spectral energy distribution (SED) 
of NGC 1068 as a template. For this SED, the optical (750nm) 
to 100 keV luminosity is a factor 11.8 greater than that in the 
20keV-100keV range (Pier et al. 1994). For NGC 2992, the lu- 
minosity in this latter range is (1 - 4) • lO"'-' erg s"' (Beckmann 
et al. lSOOTl l. Thus flie total luminosity of NGC 2992 is about an 
order of magnitude greater, which yields about 10'*'' erg s"' . This 
is a factor of 3-30 times greater than the energy estimated above 
that is required to generate the outflow. We therefore conclude 
that it is the AGN, rather than the starburst, that has driven the 
outflow. 



5.2. Velocities along the major axis 

Longslit spectra were taken with ISAAC along the major axis 
(position angle PA = 34°) at 2. 1 //m. The resulting velocity curve 
of H2 l-OS(l) along the slit shows a steep increase to the south- 
west (decline to NE) up to a radius of 1'.'5 away from the centre 
followed by a sharp decline within about 0'.'5 and subsequently 
a shallower incline. To the northeast the velocity curve behaves 
similar This is also observed in the SINFONI data (Fig.|5]l, and 
is in agreement with published data of Marquez et al. ( 119981 ) for 
a position angle of 30° and Veilleux et al. (|2001), too. 
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Fig. 5. Radial velocities (top) and corresponding velocity disper- 
sion (bottom) from longslit spectra taken with ISAAC along the 
major axis of NGC 2992 (crosses). Overplotted is a radial veloc- 
ity curve drawn from SINFONI data also along the major axis. 
1" correspond to 150pc, northeast is to the left, southwest to the and above 1" are also given (dashed line), 
right. 



Fig. 6. The dependence of the position angle on the distance 
from the centre for the 2-dimensional stellar (top) and H2 1- 
0S(1) data (bottom). The mean values of the PA for radii below 



Colina et al. (119871) found asymmetric line profiles for a re- 
gion with r < 3", and symmetric profiles for r > 3". They in- 
terpreted their results as dynamically decoupled nuclear and off- 
nuclear regions. Their velocity curve of [Oiii] is in agreement 
with our velocity curve. However, due to their coarse sampling 
along the major axis they could not resolve less than 3". 

Evidence for dynamical decoupled nuclear and off-nuclear 
regions comes also from kinemetry (Krajnovic et al. 120061 1 of 
our 2-dimensional stellar velocity data. With the position angle 
and inclination as free parameters the PA changes from 45° + 1 1° 
to 22° + 5° at about 1" (Fig.|6l). This is in good agreement with 
the 22.5° of JaiTett et al. (2003) derived from the K^-band 20 
mag arcsec"^ isophot. The PA of 45° better fits to the inclination 
of the accretion disk of about 37° and 46°, derived by Gilli et al. 
( I2OOOI 1 with the assumption of both Schwarzschild and Kerr met- 
rics from disk line models and iron line profiles, than to the 34° 
of Garcia-Lorenzo et al. (2001). The PA of the 2-dimensional 
H2 l-OS(l) velocity data also changes at 1", however from 8°+4° 
to 20° +4° (Fig.|6]). That the changes in PA are derived from kine- 
matic rather than photometric data, argues that they are real and 
not just an effect of extinction. It therefore seems that the stars 
and the gas show a different behaviour in the very centre and a 
similar one outside. 



5.3. A geometric model for the inner 3 " x 3 " 

In order to investigate the different components contributing to 
the velocity field and to explain the different position angles ob- 
served in the central 3 " x 3 " of NGC 2992 we developed a 
simple geometric model. It consists of an inner and an outer 
disk, in which particles were assumed to be in circular motion 
in a gravitational potential. The radius of the inner disk is de- 
termined by the change in PA found in kinemetry of stellar and 
H2 l-OS(l) data at about 1". It is not related to the starburst nor 
does it reflect the sphere of gravitational influence of the black 
hole, which is only a few parsec in radius. The radius of the 
outer disk is cut at 5" due to our FOV of 3". In addition, a coni- 
cal outflow with constant velocity represents the figure-of-8. All 
components could be rotated by arbitrary angles around the x-, 
y-, and z-axis, which points to the observer, thus allowing for 



different orientations to each other Velocities were adjusted in 
such a way, that they reflect the observed velocity curve. Due to 
the simplicity of the model, we only can draw some qualitative 
conclusions. 

There is an inner and outer region which must have different 
orientations in space (realized in the model by different rota- 
tion angles around the x-axis) in order to explain the different 
position angles. However, the inclination angle of 70° and the 
position angles cannot be modelled exactly at the same time. 
Dependent on the start value for the ellipticity of the unrotated 
outer disk we get an inclination angle of 50° for a circular disk 
and 63° for an ellipticity of 0.7, and position angles of 13° and 
18°, respectively. For the inner, circular region we get an incli- 
nation angle of 60° and a position angle of 34°. The different in- 
clination angles might point to a dynamical decoupled inner and 
outer region as it was already suggested by Colina et al. ( 1987P . 

In order to achieve a position angle of the figure-of-8 con- 
sistent with the pubHshed -26° (Wherle et al. 1988), it has to be 
perpendicular to a plane whose inclination angle is 51°, which 
is slightly lower than the inclination of the inner disk. One can 
identify this plane with the plane of the accretion disk, which it- 
self is not modeled. The resulting position angle is -27° with the 
SB part of it pointing away from us, consistent with the reshifted 
emission there seen by Garcia-Lorenzo et al. (1200 11 1. 

The derived geometry from our model yields the values from 
the literature of 15° for the position angle of the outer region 
(photometric axis), 34° for the inner 5" (Garcia-Lorenzo et al. 
120011), and flie inclination of the accretion disk of 46° ± 7° (GilH 
et al. '2000) very well. The resulting velocity field from this ge- 
ometry is shown in Fig.[8]and is in good agreement with the mea- 
sured velocity field of H2 l-OS(l) depicted in the lower panel of 
Fig. |2] Within this geometry the hump in the velocity curve at 
1'.'5 is the signature of the cone describing the figure-of-8, which 
is superposed on the inner disk. The stellar velocity curve along 
the major axis from SINFONI and ISAAC does not show these 
humps at 1'.'5 which speaks in favour of this interpretation. Our 
model suggests that the outer disk, inner disk, and accretion disk 
are at different orientations. One possible reason for this is that 
this geometry reflects a single warped disk, or even a warped 
outflow. However, it is beyond the scope of this paper to explore 
all the options and to assess the physical reason behind the ob- 
served geometry. 
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Fig. 7. Cartoon which shows the geometry of the central region 
of NGC 2992. Those parts of the disks with solid lines point 
towards us, those with long-dashed away from us. Short-dashed 
lines indicate the position angles of 15° for the outer and 34° for 
the inner disk. 
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Fig. 8. Calculated velocity field in the central region of NGC 
2992 resulting from two inclined disks and a conical outflow 
(see subsection 5.3.). 



6. Conclusions 

We have presented a detailed analysis of near infrared K-band 
spectra complemented by N- and Q-band spectra of NGC 2992. 
K-band spectra were obtained with the adaptive near infrared in- 
tegral field spectrograph SINFONI and allow a reconstruction 
of the 2-dimensional distribution and kinematic of the stars and 
gas in the inner 3" x 3" (450pc) at an angular resolution of (X^. 
The N- and Q-band data were obtained with Spitzer and taken 
from the archive. We compared the equivalent width of Bry and 
the supernova rate, which was derived from radio data from the 
literature, to STARS evolutionary synthesis models and find ev- 
idence for a short burst of star formation (40 - 50) Myr ago. 

From our near-infrared data as well as equivalent widths and 
luminosities of fine-structure emission lines and polycyclic aro- 
matic hydrocarbon features detected in the N- and Q-band spec- 



tra we were able to estimate the nuclear star formation rate and 
quantify the contribution of the star formation to the AGN lu- 
minosity: the luminosity is dominated by the AGN activity with 
a contribution of only 10% - 20% from nuclear star formation. 
This is in agreement with the conclusion that part of the Bry 
emission cannot be attributed to star formation but is due to the 
AGN. 

Energy and timescale considerations let us conclude, that the 
starburst would have to be more recent and more intense to drive 
the figure-of-8 as well as the large scale (2kpc) outflow. On 
the other hand the estimation of the luminosity of the AGN of 
10"*^ erg s"' is about an order of magnitude higher than is neces- 
sary to drive the outflows. 

The observed velocity curve in H2 l-OS(l) of NGC 2992 can 
be described as the superposition of the galaxy rotation and a 
conical outflow with constant velocity. The different position an- 
gles of the inner and outer region can be modelled by disks with 
different orientations in space and thus being probably dynami- 
cally decoupled. However, these disks are not related to the star- 
burst. In this picture the south east cone is pointing away from 
us supporting the findings from Garcia-Lorenzo (2001). 
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